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Abstract 

We investigate the meson, baryon and anti-baryon productions in np and Kp 
collisions in terms of a constituent cascade model which includes charm flavour. The 
constituent diquark in an incident baryon has a tendency to get a large momentum 
fraction, resulting in considerable difference of subenergies for particle productions 
from quark-diquark and quark-antiquark chains in a meson-baryon collision. As a 
consequence, the leading particle effects on D meson productions are expected to 
be quite different in 7r + and tt~ beams at present accelerator energies. 

pacs number: 13.85.Ni 



Submitted to Zeitschrift fur Physik C 



1 



Experiments at CERN (WA82) and Fermilab (E769) measured charmed meson dis- 
tributions and observed a strong asymmetry between leading and non-leading charmed 
meson productions in 7r-nucleon collisions [I], 0. The leading particle contains the same 
type of quark as one of the valence quarks in the incident hadron, while the non-leading 
one does not contain projectile valence quarks. In perturbative QCD at leading order, 
the factorization theorem predicts that c and c quarks are produced with the same dis- 
tributions and then fragment independently. Thus, the perturbative QCD model can not 
expect such a large asymmetry and contradicts with the experimental data ||. Even in 
the case of next to leading order, the predicted asymmetry is small ||. Alternatively, 
the final interaction of charm quark with spectator quarks may be responsible for the 
large asymmetry. The Lund string fragmentation model implemented with this mech- 
anism, PYTHIA Monte Carlo program, reproduces the trends of increasing asymmetry 
with Feynman x, but the asymmetry is too large |J. Recently, Vogt and Brodsky took 
account of an intrinsic cc in projectile hadron state and predicted an asymmetry from 
coalescence of charm quarks with the comoving spectator quarks, in agreement with the 
data ||. The charmed hadron production was also discussed in the framework of the 
dual string model J7|. At the present stage, the asymmetry problem about the charmed 
hadron productions is in an open question. 

On the other hand, soft hadron productions have been extensively investigated in 
terms of constituent quark chain or constituent string models [8-15]. A constituent quark 
(diquark) is a quark-gluon cluster consisting of the valence quark (diquark), sea quarks 
and gluons. Constituent quarks and diquarks play important roles in soft hadron interac- 
tions in contrast with parton quarks and gluons in hard interactions. We have proposed 
a covariant quark-diquark cascade model and applied it to various types of hadron pro- 
ductions in hadron-hadron and hadron-nucleus collisions, successfully |16|. General back- 



ground, ideas and consequences of the cascade and chain type models are reviewed in ]T] 
In our model, leading and non-leading particle effects on hadron spectra are naturally ex- 
plained by the constituent cascade mechanism. The momentum fractions of constituent 
quarks (diquarks) in projectile hadrons are characterized by parameters determined from 
Regge intercepts 0, |9|. When two incident hadrons collide, they have a tendency to 
break up into fast and slow constituents. Two cascade chains are exchanged between the 
constituents in the beam and target hadrons. The leading particle effect results from the 
first cascade step of the energetic valence constituents. It may be valuable to introduce 
charm flavour into our model and investigate the strong asymmetry between leading and 
non-leading charmed meson productions. 

In this paper, we calculate the non-charmed and charmed hadron productions in up 
and Kp collisions in terms of the covariant quark-diquark cascade model. The model 
is found to be able to explain the strong asymmetry between leading and non-leading 
charmed meson productions. Furthermore, it predicts different leading particle effects in 
D meson productions between tt + and n~ beams. In Sect. 2, we formulate the model. In 
Sect. 3, the numerical analyses of non-charmed and charmed hadron spectra are given. 
Conclusion and discussion are given in Sect. 4. 



1 Formulation of the model 



We consider a covariant cascade model 0] in terms of constituent quark and diquark. We 
assume that a baryon is composed of a constituent quark and a constituent diquark, and 
a meson is composed of a constituent quark and a constituent anti-quark. Hereafter, we 
simply call them as a quark and a diquark. We construct a Monte Carlo event generator 
to satisfy the energy-momentum conservation condition and apply it to hadron-hadron 
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collisions. We consider an inclusive reaction A + B —> C + X . In the centre of mass 
system of A and B, the light-like variables of A and B are defined as follows: 



P i Pcm B ^ T Pcm 



'so V s o 
where y^so is the centre of mass energy of the incident hadrons A and B. 



1.1 Dissociation of incident hadrons and formation of the chains 

When the collision between A and B occurs, the incident hadrons break up into con- 
stituents with a probability (1 — P g i); otherwise they emit wee gluons with P g i fol- 
lowed by a quark-antiquark pair creation. We assume four reaction types: a) non- 
diffractive dissociation, b),c) single-diffractive dissociation and d) double-diffractive dis- 
sociation types as shown in Fig. [I]. The probabilities of these types to occur are 
(1 — Pgi) 2 ,P g i(l — P g i),Pgi(l — Pgi) and PjJ , resp ect ively. Here we denote the quark- 
antiquark pair emitted from A (B) via the wee gluons as Ma ( Mb)- The probabilities 
of Ma ( Mb) to be uu, dd, ss and cc are denoted as P uU , P^ Pss and P cg , respectively. 
From the isospin invariance, we put P uU = P d ^. At very high energies, it is expected that 
In cascade chains are generated in accordance with n Pomeron exchange. However we 
neglect here the processes with n — 2, since the incident energy to be considered is not 
so high. 

In the non-diffractive dissociation, both incident hadrons A and B break up into 
two constituents. Two cascade chains are exchanged between A and B. Each cascade 
chain is developed from the two constituents, one of which is in the beam and the other 
in the target hadron. In the case of meson-baryon collisions, they are the quark-antiquark 
and the quark-diquark cascade chains as shown in Fig. la). As shown in Ref. ||, the 
distribution functions of the constituents in the incident hadron A composed of a and a' 
are described as 

H a/A (z) = H al/A (l -z)= ) \ . (2) 

Then the light-like fractions of a and a' are x\ = x$ + z, xt = Xq_R, = Xq + — 

respectively. Hereafter R denotes the uniform distribution in the 

interval < R < 1. 

In the single diffractive dissociation, we assume that the diffractively dissociated 
hadron, say B, breaks up into two constituents and the other incident hadron emits 
a wee quark-antiquark pair MA^qq) and two cascade chains are exchanged between Ma 
and B. The single-diffractive dissociation mechanism results in a 
diffractive peak of spectrum of the same kind of the incident particle A at x ~ 1. We 
assume the distribution function 

H Ma/a {z) = z^-\\ - zf^/BiPgt,^), (3) 

with the light- like variables of Ma defined as, 

xX A = x£ + z, x*? A = xtR- (4) 

Then the leading particle A and the diffractively dissociated hadron B have the following 
momentum fractions: 

xX = x 4 + (l -z), x^ = m A /(xi s ), 
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x_ — Xq_ (x_ Xq_ (1 -R)), x + — x _j_, (5) 

where the mass shell condition is considered and transverse momentum of the leading 
particle is neglected. In the centre of mass system of Ma and B, we define the light-like 
variables of these hadrons and fix the light-like fractions of the projectile constituents 
from Eq.(0) as in the non-diffractive dissociation. The single-diffractive dissociation of A 
can be treated similarly, exchanging the role of A and B. 

In the double-diffractive dissociation, both the beam A and target B emit wee mesons 
Ma and Mb and then break into two constituents. The momentum fractions of Ma and 
Mb are determined by Eq. (^) in the same manner as in the case of single-diffractive 
dissociation. We have 

„M A — „A M A _ A d 

<Aj _|_ Ju Q_J_ & 5 «X q J. 



x 



^ = X B_ Z ^ X M B = X B +R ,_ (6) 



The momentum fractions of diffractively dissociated 
hadrons A and B are 



x B = xl(l-z'), x*=x* + (l-R'). (7) 

Two cascade chains are exchanged between Ma and B and between A and Mb as shown 
in Fig. Id. Also the momentum fraction of the constituents are fixed by Eq. (|2|) in the 
centre of mass system of Ma and B (A and Mb), as in the non-diffractive dissociation 
case. 

The dynamical parameters /3's in Eq. (0), which determine the momentum sharings 
of the constituents in the incident hadrons, are related to the intercepts of the Regge 
trajectories as 0. 

At = Ai = i - <v- w (o), A, = l - cty(o), 



A = 1 - aj/</>(0). 

From previous analyses |15, 16 1, we determine the values for diquarks as 



%] = 7[ ] ( A + AO > /%} = to ( A + A) > ( 9 ) 

where [ij] and {ij} denote the flavour antisymmetric and symmetric diquarks, respectively. 

The parameters j\ i and 7{} are related to the proton structure and determined from 
the ratio tt + /tt~ in proton fragmentation. The parameters (5 g i and Ad i n Eq. (ff) and P g i 
are chosen to reproduce the diffractive peaks of 7r ± X at x ~ 1 and np — > at 

x ~ — 1. The ratio P sS / P u u is fixed by the ratio n + /K + in the proton fragmentation 
region. The probability P c5 is determined from the cross section in np collisions. 



1.2 Hadron productions from the cascade-chains 

Secondary hadrons are produced from the cascade chains exchanged between the two 
hadrons: (a)A and B, (b) A and M B , (c) Ma and 5, (d) A and M B , and M4 and 
B. We consider lower lying hadrons: pseudoscalar(PS'), vector(V), tensor(T) mesons, 
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octet(O) and decuplet(P) baryons composed of u, d, and s flavours and the correspondings 
with charm flavour. We assume the corresponding meson production probabilities in the 
cascade processes as Pps, Pv an d Pt (Pt = 1 — Pps — Pv) an d the corresponding baryon 
productions as Po and Pd (Pd = 1 — Po)- Octet and decuplet baryons are described as 

|8 >= cos %[</</'] > + sin %{</</'} >, 



|10>=|g{gY}>. (10) 
We take account of the following cascade processes: 



q - M(qq') + q', B(q[q' q"}) + [q>q»] 
B(q{q'q"}) + {q'q"}, 



[q'q"\ -> B(q[q'q"])+q, M(qq') + 
M(qq') + {qq"}, 



{q'q"} -> B(q{q'q"})+q, M(qq') + [qq% 

M(qq') + {qq"}, (11) 

where M(qq") is also produced with the same probability as that of M(qq') from a di- 
quark, and meson production probabilities from q, [q'q"] and {q'q"} are 1 — e, rj\ i and 
respectively. We assume equal probabilities of symmetric and anti-symmetric diquak 
productions from a diquark for q ^ q". The pair creation probabilities P uU , P^i, P s g and 
P c5 are also assumed in the cascade process. The probabilities of [qq'\ [qq'\ , {qq'}{qq'} and 
{qq}{qq} pair creations from a quark are chosen as P qq P q i q <, P q qP q iqi and P qq - 2 , respectively. 
Then the emission probabilities of individual cascade processes are determined from the 
above probabilities. For examples, the probabilities to produce tt + , p + , aj, A ++ , p, 
A + , S+ + and from a u quark are as follows: 

(1 - e)P dd P PS , (1 - e)P di P Vl (1 - e)P dd -P T ,..., 
e PJ, e P ua P dd (1 + I P sin 2 6/( | P Q sin 2 9 + § P D ) ), 
e P ua P dd - | P D / (| P sin 2 6 + § P D ) 



' J * * * 5 



e P cc - 2 § P sin 2 ^/ (| P sin 2 £ + | P D ), e Pj \ P D / (§ P G sin 2 5 + | P D ), 

where the factors 1/3 and 1/2 are flavour SU(4) factors. 

We redefine the light-like fractions of the incident constituents in the rest frame of the 
cascade chain. The momentum sharing of the cascade process q + q' — > M(qq") + q" + q' 
from a q with x± and q' with x± takes place as follows [0, [l^|: First, using the emission 
function 

F q „ q {z) = z^-\l - zf^'^/B^ (3 q + (3 q „), (12) 

we fix the lightlike fractions of q" and Mas 4 = x q + z and x* f = x\ — x% , respectively 
and put xl = xL. Second, the transverse momentum of M is determined from the 
probability function 



GW) = ^ex P (-^Lp2) (13) 



in pip space. Then, from the onshell condition, x_ is fixed as x_ = {m M + p T )/xt;s', 
where y/s' is the subenergy of the cascade chain. The transverse momentum of q" is 
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Pt = Pt ~Pt- The lightlike fraction of q' is decreased to = scl — ar?f . If the energy of 
q' is enough to create another hadron, the cascade such as q" + q' — > q" + q'" + M(q"'q') 
takes place in the opposite side. 

The production probabilities Pps and Py are fixed from ir, p, ci2 spectra. The prob- 
abilities e and Po are chosen so as to reproduce baryon productions in meson fragmenta- 
tions. The parameters n 1 and T]i\ are related to the spectra of the baryons which contain 
no projectile diquarks in the baryon beam fragmentation region. The weight of antisym- 
metric diquark in the incident baryon cos 2 9 is chosen to get an appropriate n + /n~ ratio 
in the proton fragmentation region. By introducing the dynamical parameter 7 in Eq. 
(|T2D, we can explain the ^-distributions of various kind of hadrons by using the same 
parameters /3's in Eq. @. The parameter C in Eq. (|l^) is fixed from the experimental 
data on p\ distributions. These parameters are fixed so as to get an overall fitting to 
various hadron data with different incidents at different energies. 

1.3 Final step 

In the final step, the constituents from both sides terminate and recombine into one or 
two hadrons according to the processes: q + q' — » M(qq'), q+ [q'q"] — > B(q[q'q"]), {qq 1 } + 
[q"q"'} -> B(q*{qq'}) + B(q T [q"q"']) 

and so on. The momenta of the recombined hadrons are the sum of those of the final 
constituents and are offshell. The energy excess of the recombined hadron on i-th cascade 
chain is 

AE {i) = E' {i) - E {i) , 

E ( i ) = y Jptf + m? (i) , i = l,2, (14) 

where E't* and p^ are the sum of the energies and three momenta of the final constituents 
of i-th chain. There are different methods to put these particles on mass shell: (1) Mul- 
tiplying the three momenta of all produced hadrons by a factor a so that the summation 
E"=i \J ( a Pi) 2 + m l would be e( l ual to v^o- ( 2 ) Whenm 2 2 = (E'^ + E'^f -(p {l) +p {2) ) 2 > 
(j7i(i) + m( 2 )) 2 , we can use two body decay m 12 — > + m( 2 ). Here we use the second 
method. In order to minimize the effect of this procedure, however, we choose a pair of 
hadrons j and k which has the smallest value of m 2 - k — (rrij + mk) 2 = (Ej + Ek + AEq) + 
A£'(2)) 2 — (Pm +P(2)) 2 ~ ( m j + m k) 2 > among produced hadrons in the collision between 
A and B. In the two body decay, the directions of three momenta in j + k system are 
conserved. The energies and momenta of all directly produced hadrons are fixed. Then 
directly produced resonances decay into stable particles. For simplicity, their decays are 
assumed to be isotropic in the rest frame of the resonance. 



2 Comparison with data 

The newly introduced parameters concerning the charm flavour are the probability of 
cc pair creation P c5 and the parameter f3 c fixing the momentum sharing of the charmed 
constituent in the cascade processes. From the data on tt~p — > D^X cross section |19| , 



we choose the parameters as P c5 = 0.00016 and P sS = 0.09984 . The dynamical parameter 
(3 C is related to the Regge trajectory of J/ip. Adopting the value of a' = 0.9, we take the 
value as (3 C = 8.0. 

We set other parameters related to the u, d and s flavours to be nearly equal to the 
values used in the previous analyses |[5|, |i~6"l . Namely, the pair creation probabilities of 
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uu and dd are P uU = P d d = 0.45. Meson and baryon production ratios are chosen as 
P PS : P v ■ Pt = 0.4 : 0.4 : 0.2 and P Q : P D = 0.33 : 0.67, respectively. From p - uo and 
trajectories, (3 's are chosen as /3 U = Pd = 0.5 and (3 S = 1.0, respectively. The dynamical 
parameters in Eq.@ and in Eq.([[|) are chosen as 7[ ] = 1.5, 7{} = 2.0 and 7 = 1.75 
. The parameters P g i, (3 g i and (3id are fixed from the diffractive peak of the projectile 
hadrons at \x\ ~ 1. The probability of the incident hadron to emit a wee meson is set 
to be P g i = 0.15 and the parameters (3 g i and flu are chosen as fl g \ = 0.1 and flu = 3.0. 
(Changing the parameters (3 g i and flu do not alter the shape of the spectra so much.) 
The mixing angle of symmetric and anti-symmetric diquarks in octet baryons is put as 
cos 2 9 = 0.5. The meson emission probabilities from quark and diquark are chosen as 
1 — e = 1 — 0.07 and r/[ ] — V{} — 0.25, respectively. The p\ distribution parameter in 
Eq.flTS]) is chosen as C — 1.2 in GeV unit. 



2.1 Non-charmed hadron spectra 

In Fig. ^, the p\ distribution of 7r~ in 7r + p collision is compared with the experimental 
data at pi = 250 GeV/c |20|| . Method (1) suppresses the distribution at p\ > 5 (GeV/c) 2 
as compared with method(2). For x distributions and asymmetries, the difference between 
methods (1) and (2) is very small and in the following we show only the results of method 
(2). In Fig. |3|, we show the results of 7r ± and distributions in n + p collisions at pl = 3.7 
GeV/c [^T|] and 100 GeV/c P2"| . We note that the four momentum is conserved by our 
covariant formalism. The s-dependence as well as the x behaviour of the spectra are well 
reproduced. The kinematical constraints on particle productions reduce the spectra of 
7r~ and K~ at pl = 3.7 GeV/c relatively as compared with those at pi = 100 GeV/c. 
The enhancement of the 7r + spectrum at x = 1 is due to the diffractive target dissociation 
type mechanism. The leading particles have harder x distributions than the non-leading 
particles. From the choice of the dynamical parameters /?'s in Eqs.(||) and (|j), the incident 
u and d valence quarks in proton have softer distributions than the u and d valence quarks 
in 7r + . Therefore the spectra of leading mesons in the proton fragmentation are softer 
than those in the ir + fragmentation. Even the non-leading K~ spectrum has very similar 
tendency, reflecting the initial configuration of the constituents of the chain. In Fig. ^, 
we show the result of proton distribution and compare with the data El. The peak of 



the proton distributions at x ~ — 1 is also attributed to the single-diffractive dissociation 
type mechanism. 

The x distributions of p° and in n + p and K* + and in K + p collisions at pi = 250 
GeV/c [§4| are shown in Figs. |5] a and b, respectively. The peaks of spectra at x ~ 



are steeper than those of vector mesons in ir + p and K + p collisions, reflecting resonance 
effects on PS mesons. The spectrum of in K + p interaction at x > decreases more 
slowly with x than that of K® in ir + p interaction. Here we see the effect of the constituent 
distributions i.e. the s in K + is energetic as compared with the constituents u and d in 
7r + . The agreement with the data is fairly good. 



2.2 Charmed hadron spectra 

Next we show the result of charmed meson productions. In Fig. |6], the leading (Did) and 
non-leading ( D n i) D meson productions with respect to the incident ir~ in ir~p collisions 



at pl = 360 GeV/c are compared with the data [ 1~9| . Here we mean and D n i as 



D ld = D- + (D*~ -> D°) and D nl = D + + (D* + D°), respectively. From these data 
we fixed the value of P c£ . 
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We calculate the asymmetry A which is defined as 



D+ + d*+ - D~ - D*- 

A n ± = ± (15) 

D+ + D*+ + D- + D*- v ; 

for n ± p collisions and compare with the data |]I|, |2| in Fig. |^. The experiment WA82 
measured the asymmetry by 340 GeV/c 7r~ beams and E769 used a mixed beam (75%7r~ 
and 25%7r + ) at 250 GeV/c. Fig. |7|b) shows our calculated p\ dependence of A n ± in the 
region 0.1 < x < 0.7 with the data of E769. Our model predicts rather large positive 
values for tt^ 1 beams with A n + > A n ~. The discrepancies between the prediction and 
the experimental data may be attributed to the neglect of multi-chain type mechanism 
and/or inadequate p\ distribution. 

In 7i~p interaction, the leading particle D~ is produced on the chain between the va- 
lence d quark in the beam and the valence diquark in the target. While in n + p interaction, 
the leading particle D + is produced on the chain between the valence d quark in the 7r + 
beam and the valence quark in the target. From Eqs.(4) and (5), the proton target has a 
tendency to break into an energetic valence diquark and a wee valence quark. The former 
chain developed from d quark and diquark in n~p collision is more energetic than the 
latter chain. The leading D~ meson produced in the first cascade process is energetic. 
Since the most energetic non-leading D + is produced after D~ production, D + is less 
energetic in our model. Therefore the resulting asymmetry is large at x ~ 1. In the case 
in which a baryon is produced from the valence diquark in the first cascade process, the 
momentum of the valence d quark is reduced. Furthermore, when the valence d quark 
has a small momentum, the total momentum is shifted to the diquark side. Then the 
leading spectrum of D~ on the quark-diquark chain tends to have a small momentum 
and the asymmetry is reduced at x < 0.7. On the other hand, the d and the quark chain 
in the n + p interaction is less energetic than the quark-diquark chain. So the leading D + 
production with a large momentum in the first cascade process is suppressed resulting a 
small asymmetry A n + at x ~ 1. The momentum reduction of the d quark is small in the 
case when a meson is produced in the first cascade step from the valence quark. Therefore 
the asymmetry A n + is not so suppressed as compared with A n - at x < 0.7. 

In the diffract ive dissociation type mechanisms, the chain belongs almost to either the 
beam or the target fragmentation and particle productions have no ir^ beam dependencies. 
With the four types of interaction mechanisms the resulting asymmetry in ir~p collision 
increases with x and has a dip structure in the vicinity of x ~ 0.5. The asymmetry in 
7i + p has a peak at x ~ 0.6 and smaller than the one in 7r~ beam at 0.8 < x. 

In Fig. |8] the A and A c distributions computed in our model are shown with the A 
spectrum P5[. The resulting spectrum of A c at x < has a steeper peak than that of 



A. This is attributed to the emission function with the large value of j3 c . Leading baryon 
(which contains the same kind of diquark as the incident target) spectra have a peak at 
x < due to the valence diquark effect. It is also noted that the non leading baryon 
spectra at x < are also affected by the valence diquark in the target. The shapes of the 
non-leading baryon spectra are softer than those of leading baryons, but are harder than 
those of antibaryons and mesons in p fragmentation. This is because of the effect of the 
processes qq — » M + qq'. In Fig. |9] we show our calculations of the A and A c distributions 
and compare with the measured A data ||23|| . The A spectrum in K + beam fragmentation 
shows a strong leading particle effect as compared with that in 7r + fragmentation. This 
is because the valence s quark in the K + beam is harder than the valence d quark in ir + 
beam. 
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3 Conclusion and discussion 



The covariant quark- diquark cascade model with charm flavour has been investigated. 
The model can successfully explain the non-charmed hadron spectra. The mass effects 
on the spectra are well explained by the flavour SU(3) breaking parameters Pfj and the 
momentum sharing parameters f3j. As a result of the different hardness of the valence 
constituents in the incident hadrons, the leading particle effects on hadron productions are 
vastly different from each other. The order of the hardness of hadron spectra characterized 
in view of the constituent quark and diquark cascade is as follows: in decreasing order are 
p — > B, M — > M,p — > M, M — > B(B) and p — > B. The detailed features of the spectra 
are described by the flavour dependencies of Pfj and j3j, and by reaction mechanisms in 
Fig. |. 

We have investigated charmed hadron spectra. The dynamical parameter f3 c deter- 
mined by Eq.(§) and P c5 well explain charmed hadron productions. It is noted that the 
dual-Regge ansatz for constituent quark distributions [§] appear to work even for charmed 
quark. The momentum sharing functions in cascade processes and the beam dissociation 
functions are characterized by the same dynamical parameters determined from Regge 
intercepts. Energy conservation condition with the constant values of Pn well explains 
various hadron production ratios in the energy range 3.7 GeV/c < pi < 360 GeV/c. 

In our model, the asymmetry between leading and non-leading meson productions 
depends on the incident beams. The asymmetries A n ± in Fig. |7| reflect the distributions of 
d in 7r + and d in ir~ beam. Although the distributions of the constituents in the incident 
7r =1= are the same, we have different asymmetries between ir + and 7r~ beams. In the single 
and double-diffractive dissociation, particle productions in the beam fragmentation are 
separated from those in the target fragmentation region. However in the non-diffractive 
dissociation, the particle productions are mixed up with each other. We have different 
cascade chain properties between the quark-antiquark and the quark-diquark chain and 
have different leading particle effects between 7r + and ir~ beams. Therefore the weight of 
the non-diffractive dissociation mechanism (1 — P g i) 2 is also related to the the asymmetries 
A with respect to 7r =l= beams. 

Although it is interesting to apply the model to hadron productions containing bot- 
tom quark, their cross sections in soft interactions are too small to discuss seriously at the 
present stage. So we disregard the possibility of the hadrons containing bottom quark as 
a valence constituent. 
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Figure 1: The reaction mechanisms in AB collisions Non-diffractive dissociation type, b,c 
Single-diffractive dissociation type and d double-diffractive dissociation type mechanisms. 



Figure 2: Calculated values for p\ distributions of n in n + p collisions at p L = 250 GeV/c. 
Experimental data are taken from [20] . 



Figure 3: Invariant x distributions of n ± and K ± in n + p collisions at a p L = 3.7 GeV/c 
and b p L = 100 GeV/c. Experimental data are taken from [21] and [22]. 



Figure 4: Calculated values for x distributions of p in ir + p collisions at pi = 250 GeV/c. 
Experimental data are taken from [23] . 



Figure 5: Feynman x distributions of a p° and K$ in ir + p and b Kg and K* + in K + p 
collisions at pl = 250 GeV/c. Experimental data are taken from [24]. 



Figure 6: Calculated values for x distributions of a D\ d = D~ + (D*~ — > D°) and b 
D nt = D + + {D* + — > D°) in n~p collisions at p L = 360 GeV/c. Experimental data are 
taken from [19]. 



Figure 7: Asymmetries of D ± in ^p collisions: a vs. x and b vs. p\ at 0.1 < x < 0.7. 
Experimental data are taken from [1] and [2]. 



Figure 8: Feynman x distributions of A and A c in a ir + p and b K + p collisions at pi = 250 
GeV/c. Experimental data are taken from [23]. 



Figure 9: Feynman x distributions of A and A c in a ir + p and b K + p collisions at pl = 250 
GeV/c. Experimental data are taken from [23]. 
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Fig. 3 
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